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I ■ Abstract 

^ \ Frequency sum rules are derived in extended quantum systems of 

non relativistic fermions from a minimal set of assumptions on dynam- 
ics in infinite volume, for ground and thermal states invariant under 
space translations or a lattice subgroup. 
^ ■ For the jellium Coulomb model, they imply the one point result for 

J;^ . the plasmon energy spectrum in the zero momentum limit. 

\^ I In general, the density waves energy spectrum is shown to converge, 

^^ ' in the limit of large wavelenght, to a point measure at zero frequency, 

T^ljj- ■ for any number of fermion fields and potentials V^"^^' with integrable 

O I second derivatives. 

For low momentum, < uj^{k) >~ k"^ for potentials with r^ didjV inte- 
grable, < uJ^{k) >~ k°'~ '^'^ for potentials decaying at infinity as 1/r", 
d — 2<a<d, d the space dimensions. 
K> ' For one component models with short range interactions, the fourth 

C^ • momentum of the frequency is expressed, at lowest order in k, purely 

in terms of the three point correlation function of the density. 



Math. Subj. Class.: 82B10, 82B21 
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1. Introduction 

The physics of charge density waves in extended systems was under- 
stood in ref.jl] in terms of classical charge configurations giving rise 
to slowly varying electric fields, resulting in plasma oscillations with 
frequency ujp = e^p/m for infinitely extended systems of particles of 
mass m, charge e and density p. 

In solid state physics such an analysis applies [2j on the basis of 
suitable simplifications, in particular oi a, random phase approximation 
[3] in the analysis of charge density correlation functions. Plasma 
oscillations correspond to a single point frequency spectrum, uj = Up, 
for the charge density correlation function at zero momentum, and it 
is not clear whether neglected terms may spoil such pure point result, 
even in the case of only one kind of charged particles in a neutralizing 
uniform background (jellium model). 

The control of the energy spectrum associated to density waves also 
plays a crucial role in the the theory of quantum liquids, in particu- 
lar in the discussion of superfluidity, [1] [S] [S]; very similar problems 
appear in the discussion of the current commutators at low momen- 
tum which are at the basis of the Kubo and Landau approaches to 
superconductivity [7] [8]. 

In ref. [9] the relevant spectrum in the jellium Coulomb model was 
derived, and shown to consist only of the plasma frequency, through 
an analysis of Galilei transformations, which are spontanously broken 
in the jellium model, as in any nonrelativistic system at non-zero 
density. 

The derivation employed a generalized Goldstone theorem, which 
gives the energy spectrum associated to spontaneosly broken symme- 
tries in terms of commutators between generators of symmetries and 
order parameters. Such commutators are time independent in the 
ordinary Goldstone case, but not in the presence of sufficiently long 
range interactions, the Coulomb potential in the jellium case. Similar 
results have also been obtained in refs. [lOj through an analysis of 
operators describing long wave charge fluctuations in the jellium and 
other models. 

The Galilei analysis is complicated by the lack of invariance un- 
der space translations of the density of the Galilei generator and by 
short distance problems associated to the singularity of the Coulomb 
potential at the origin. The analysis can in fact be simplified by focus- 
ing on charge density commutators, a tecnique discussed e.g. in refs. 
|llj . [12j . [13j . On one side, such commutators provide an alternative 
derivation of the spectrum associated to Galilei transformations, on 
the other they have convenient positivity properties and simple evo- 



lution equations. Moreover, as we shall see, the charge commutators 
analysis applies to a large class of systems and allows for the derivation 
of exact relations for (density waves) energy spectra, also at non-zero 
momentum. 

The aim of this paper is to present a self contained derivation of 
the plasma spectrum at zero momentum in the jellum model, based on 
the analysis of charge commutators, to extend the analysis to the case 
of spontaneous breaking of translation to a lattice subgroup (which is 
relevant in presence of Wigner crystallization [14j, ^15]) and to derive, 
by the same methods, general relations for the frequency spectrum at 
low momenta of density waves in non relativistic fermion systems with 
short and long range potentials. 

The use of sum rules to obtain information on the density waves 
frequency spectrum seems to have been prevented by results on the 
divergence of momenta of order > 3 [16]. However, the divergence of 
the third momentum only appears for singular (delta or hard core) 
potentials and the divergence of the fifth momentum in the Coulomb 
case depends on the singularity at the origin of the potential. More 
generally, the (perturbative) results of [TB] indicate that no momen- 
tum diverges for regular (C°°) potentials, a result which also follows 
from our analysis. In the following, the Coulomb potential will be 
regularized at the origin and the limit A; — > of the energy spectral 
measure will be shown to be independent of the regularization. 

We remark that the result on the plasma spectrum at zero momen- 
tum does not follow from the sum rules for the frequency momenta 
up to the third, nor by the use of the "perfect screening sum rule" , 
concerning (cj^^) [6], unless a one point approximation is assumed for 
the frequency spectrum; however, the reduction of the specturm to a 
single point is the key result, which involves the fifth momentum and 
holds, in our analysis, only in the limit fc — > 0. 

Besides the analysis of the plasma spectrum and of the fifth fre- 
quency momentum, we reobtain the third momentum sum rules of 
refs.[Tl] [T7j [12], clarifying and symplifying their derivation, extend- 
ing them to states with discrete translation symmetry, and deriving 
exact consequences on energy spectra, in particular in the case of long 
range potentials with a faster decay with respect to the Coulomb po- 
tential. Our results do not require the explicit construction of ground 
and thermal states and may also shed light on the problems and al- 
ternatives which appear in their analysis [TH], [19j . |2Uj . |21j . 

We consider infinite systems of nonrelativistic fermions described 
by the canonical anticommutation relations (CAR) algebra, in d space 
dimensions, in particular d = 2,3. The time evolution is assumed to 
be given, through equal time commutators, by a free Hamiltonian and 



interaction potentials y(|a;|). We consider states Q^, with nonzero 
mean particle density (p) , invariant under space translations or under 
a lattice subgroup of them, time evolution and parity, satisfying either 
the spectral condition, i.e. positivity of the energy in the resulting 
representation of the CAR algebra, or the KMS condition at inverse 
temperature (3. 

We denote by .F(/) or by / the Fourier transform of /, hy f * g 
the convolution of / and g, by A(f) the smearing of the distribution 
^4(3;) with /; unless specified differently, test functions belong to the 
Schwartz space <S of infinitely differentiable functions of fast decrease. 
For simplicity, the Planck constant (divided by 2ir) is omitted and 
the variable io is used for the spectrum of the Hamiltonian; sum over 
repeated indexes is implicit and dx will denote d x. Our results are 
the following: 

i) one point plasma spectrum for jellium at zero momentum: 

In the jellium Coulomb model, i.e. for non-relativistic fermions in 
three space dimensions with interaction V{x) = e^ /^'k\x\ * 7(|x|), 
7 G 5 real, J ^{x) dx = 1, we consider the expectations 



d,jif{u:)= j{p{f)dE{u)p{f)) , 

with f{x) = e-x.'pikxaR{x) (see eq.©, dEiuj) the spectral measure of 
the Hamiltonian, ( ) the expectation on a thermal or ground state, 
invariant under space translations or a lattice subgroup of them and 
parity. The normalized positive measures 

dvfibj) = Nf uj [dpf{ijj) — dp,f[—uj)) . 

converge, for ur — > 1 and then A; — s- 0, to 

1/2 {6{LO-ujp) + 5{u: + u:p)) , u^l = e^{p)/m (1) 

independently of the regularization 7. The derivation relies purely on 
Newton equations. Gauss' law at large distances and invariance under 
(a lattice subgroup of) translations. 

ii) Goldstone spectrum for potentials with integrable second 
derivatives: 

For any number of fermion fields and all regular potentials V^^^> with 
didjV^ '{x) integrable, the above measures dvf, defined by the total 
density p = YliPii converge to 5(a;), correponding to a Goldstone 
spectrum (also associated to the spontaneous breaking of the Galilei 
charges Y.i I dx pi{x)xk )• 

iii) dependence of the spectrum on the decay properties of 
the potential: 



In the same cases as in ii) and with the same notation, for R ^ oo, 
i.e / -^ expifcx, 



cvUufiuj) ^c\k\^{l + o{\k\)) 
if x^didjV^x) is integrable; 

I u)'^dvf{uj) <c\kf 



if \x\f^didjV{x), < /? < 2, is integrable, with an equahty, up to o{\k), 
for one fermion field and V{x) ~ l/|x|", a = (3 + d — 2. In particular, 
< uj^{k) >= 0{\k\) for jellium Coulomb systems (i.e. V ~ l/|a;|) in 
two space dimensions. 

iv) fourth momentum of the spectrum at the order k'^: 
For one fermion field, regular potentials of fast decrease and states 
invariant under space translation, rotations and parity the fourth mo- 
mentum of diyf{Lo) converges, for / — > expikx, to 

k^/{p)m^ j dxdyW{x,y){p{^)p{x)p{y)) + 0{k^) , (2) 

with W{x,y) computed in Sect. 6 in terms of the potential. In partic- 
ular, < Lo'^{k) > is not of the same order as < uP'{k) >^, as it would 
follow from a single quasi-particle interpretation. 

In Sect. 2 the mathematical framework is specified, in terms of 
correlation functions for infinite systems, with time derivatives given 
by appropriate commutators with local hamiltonians. Energy spectra 
are expressed in terms of time derivatives of commutators; conver- 
gence to delta functions of the spectral measures in the limit of zero 
momentum follows from relations between frequency momenta. 

In Sect. 3 equations of motion, commutators and low momentum 
expansions are discussed in general, together with their implications 
on energy-momentum spectra for short range interactions. 

In Sect. 4 the expression for < uP'{k) > is shown to result, for long 
range potentials, in the above dispersion relations. 

In Sect. 5 the fourth momentum of the frequency is computed in 
the jellium model in the limit fc — > 0, implying the plasma frequency 
result. 

In Sect. 6 the first term in the k expansion of < uj^{k) > is 
expressed in terms of the three point function of the density for regular 
short range potentials. 



2. Commutators and energy spectra 

We consider infinite systems, described by fermion fields il)*{f),ilji{f)^ 
generating an ACR algebra [22], with 

[i;*if),ijj{g)U = 6ijJf{x)g{x)dx , [Mf),i^ji9)]+ = , (3) 

dx staying, here and in the following, for d'^x in d space dimensions. 
The addition of spin indexes leads to minor changes, which will be 
indicated when relevant. The results of Sects. 3 and 4 hold for any 
number of fermion fields, those of Sects. 5 and 6 in the case of only 
one fermion field. 

We consider representations defined (through the GNS construc- 
tion [22j , \23\ ) by states invariant under (a lattice subgroup of) space 
translations, with a finite number of particles when restricted to finite 
regions, i.e. defining locally Fock representations. This allows for the 
use of variables in the weak closure of the ACR algebra in the Fock 
space, e.g. bounded functions of the operators p{f), (density opera- 
tors integrated with regular functions), and in fact we will work, in 
the spirit of Wightman theory |24j . with unbounded field operators 
likep(/). 

In the presence of interactions, the construction of the dynamics 
of such systems in terms of automorphisms of the ACR algebra is not 
completely under control; in the case of lattice spin systems [25] , inte- 
grability of the interaction implies norm convergence of finite volume 
dynamics and stability of the quasi-local algebra, while the results 
of [26] and [9] imply the necessity of weaker convergence and larger 
algebras in the case of Coulomb interactions. 

For our purposes, it is enough that dynamics exists as a group of 
automorphisms of an algebra containing the ACR algebra and that 
expectation values of time derivatives are given, at zero time and for 
suitable variables, by limits of commutators with hamiltonians asso- 
ciated to a sequence of bounded regions invading the space. We end 
therefore with the following assumptions, in the spirit of Wightman 
theory: 

A) The correlation functions at all times of the fermion fields and of 
their Wick ordered polynomials are distributions in the space vari- 
ables, continuous in time as distributions, invariant under time trans- 
lations and space translations, or under a lattice subgroup of the latter, 
satisfying Wightman positivity and therefore defining operator valued 
distributions on the invariant Wightman domain, with fundamental 
vector ^0 • The corresponding state (invariant under time and space, 
or lattice, translations) will be denoted by fi, expectations on J7 by ( ) 
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In the space translation invariant case, the unitary groups implement- 
ing space and time translations are then strongly continuous; their 
generators will be denoted by P and H, and the corresponding joint 
spectral measure by dE{k, u). For states invariant under a lattice sub- 
group of space translations, dE{k, u) will denote the spectral measure 
associated to time and lattice translations, with k in the fundamental 
cell of the reciprocal lattice; dE{uj) will denote the spectral measure 
oiH. 

B) i7 will be assumed to be a ground or a thermal state, i.e. to satisfy 
either positivity of the energy i7 > 0, or the KMS conditions 

^p{AB{t))=np{B{t-iP)A) (4) 

with B{t — i[3) defined by analytic continuation of the correlation 
functions. The correlation functions at equal time will be assumed 
to be bounded, after smearing, with respect to translations of space 
variables, to satisfy the cluster property and to be invariant under 
parity. Invariance under rotations will be assumed only to simplify 
the results, in the space translation invariant case. 

C) The correlation functions of time derivatives, at equal times, of 
Wick polynomials <I>(x) are given by the infinite volume limit of the 
correlation functions of their (multiple) commutators with the local 
hamiltonians 

(-i)"d"M"$ = [ii"... [if, $]...] = 

= lim hm [Hl,,rA...[Hl^,R^,^...]] (5) 

Li—foo Hj—KX) 

with 

= l/2m / dx aji{x) diil)*{x) diil){x) + fi dx aji{x)ip*{x)il){x) + 

+ 1/2 dxdy aR{x)aR{y)ij*{y)ij*{x)VL{x -y)ij{x)i;{y) , (6) 
with Vl{x) = aL{x) 7 * y(|a;[), 7([a:|) G S real, f "fd'^x = 1, and 

Q/j(x) = a{\x\/R) , a G 5 , a{x) > 0, a(0) = 1 , a(x) = Vx > 1 

(7) 
and similarly for more than one fermion field; the limits Rj — > cxo exist 
by locality of the ACR relations; existence, and independence of the 
order, of the limits Lj — > 00 follows from integrability of V{x) and 
boundedness of the correlation functions in the above sense. 



For gauge invariant variables, the only ones to be considered in 
the following, the second term in eq. ([6]), as well as a possible term 
pVi* (p), corresponding to an interaction with a uniform background, 
are irrelevant in the above commutators; moreover, the integrability 
condition applies to diV{x), since only derivatives of V appear, as 
a consequence of the vanishing of the commutators between gauge 
invariant variables and the integral of the charge density. 

For Coulomb interactions, in order to perform the limits Lj ^ oo, 
the truncated correlation functions will be assumed to decay, after 
smearing with test functions, as nr=i I^<t(j) ~ ^a{i+i)\~^~^ ^^^ ^-H Per- 
mutations a, for some e > 0. In general, similar decay properties 
are required in the absence of integrability of diV. The ultraviolet 
regularization is necessary, in general, for the existence of frequency 
momenta. It will be omitted in the notation, writing V for j *V. 

In the above framework, we will derive constraints on energy spec- 
tra which arise directly from the equations of motion. Goldstone the- 
orem can be regarded as one of them, and follows in fact [26j from in- 
variance of the equation of motion under a symmetry which commutes 
with space translations and is spontaneously broken, under sufficient 
locality properties of the time evolution. For long range interactions, 
the latter property may fail, and in this case the proof of Goldstone 
theorem gives a relation between an energy spectrum at zero momen- 
tum and the time dependence of the expectation value of the commu- 
tator between an order parameter A and a charge operator Q^, in the 
limit i? ^ cx) [23 [26] 

In this paper, similar information on energy spectra will be ob- 
tained from an analysis of time derivatives of commutators of the 
form {[A*{x), A{y, t)]). If A ^o is in the domain of _?/", then 

{-i)''r/dr{A*A{t))t=o = f uj'' {A* dE{uj)A) . 

The energy spectral measure of the state ^^o will be denoted by 

dfj-Aiuj) = {A* dE{uj)A); the measures 

di^Ai^) = ^ {dp.A{^) - dfj.A*{-uj)) (8) 

satisfy, under the above domain conditions, 

duAiu;) a." = {-ir+'d^+'/df'+' {[A*,A{t)])t=o , (9) 



For ground and thermal states the measures di^A are positive. If 
dfiA*i^) = dfiAi^), which will follow in our case from parity invari- 
ance, they are even and determine dfiA up to 6{uj): 
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Lemma 1 If Q is a ground or KMS state and A hermitean, dvA{^), 
defined by eq.(B^, is positive. If d^A* = d^.A, then dvA is even and 
determines d^A apart from multiples of 5{uj); in this case, if dvAn 
converge as measures, dfiAn converge as distributions, apart from S{uj) 
terms. 

Proof: For ground states, dfi{uj) has positive support, so that du is 
clearly positive and determined up to c6. For KMS states, dn{—ijj) = 
exp {—f3uj) dfi{uj), so that 

di^iio) = a; (1 - e-^^) dn{io) , 

which implies positivity of dv and uniqueness of dfj, apart from mul- 
tiples of 6{u;). Convergence in the sense of distributions of the odd 
part of dfj.A„ follows from eq.dH]) and implies, by the ground state 
or the KMS condition, convergence of dfiA„ as distributions, on test 
functions which vanish at u; = 0. D 

Energy spectra as functions of space momentum (pseudo-momen- 
tum in the case of lattice invariance) depend on commutators of the 
form {[A*{x,0)A{y,t)]); for A hermitean and of definite parity, and 
f{x) = f{—x), i.e. / real, parity (P) invariance of the dynamics and 
of implies 

{A{f,0)A{f,tr) = {PA{f,0rA{f,t)P) = {A{f,0rA{f,t)) , (10) 

i.e. dfj,A{f)*{ijj) = dnA(f)i'->j), so that Lemma 1 applies and the Fourier 
transform 

duAif) (iv) ^ T i-id/dt {[Aif, 0),Aif, t)]) ) 

is a positive even measure. For states invariant under space transla- 
tions, 

{[A{f,Or,A{f,t)]) = I {A{dE{k,io) - dE{k,-u;)) A) f{kf e^'^' , 

(11) 

so that 

di^Aif){^) = ^C^ - e~^'^) l'{AdE{k,u;)A)f{kf (12) 

for a KMS state, and similarly for a ground state. For fk,R{x) = 
e^^^ aR{x), Or as in eq.([7|), R -^ oo, 

duAj.Juj) ^ l/R" dvAUu.B)^^) ^ ^(1 - e"^") {AdE{k,u:)A) (13) 

apart from a constant, as distributions in u), assuming a polynomial 
bound in t for the L^ norm of the commutator ([^(0,0), ^(y,t)]) ^. 
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For states invariant under a lattice group of space translations, 
{[A{x + a, 0)A{y + a, t)]) is periodic in a and, for / = fk,R, apart from 
an irrelevant constant 

l/RU[A{f),A{f,t)])^ 

dx dy J^([A(x, Q),A{y + n, t)]) e^''^^'^^ e'^"" = 
AxA „ 

dx dy f{A{x, 0) {dE{k, uj) - dE{k, -u)) A{y, 0)) e^'^^y-''^ e'^* 



'AxA J 

for R — > CO, assuming integrability of the above commutator, with A 
denoting a lattice cell (the integrand being periodic both in x and y). 
Therefore, as distributions in to, 



dvAjk ijC"^) ~^ const uj{l — e 



~l3uj\ 



X 



dxdy{A{x,f))dE{k,uj)A{y,{)))e'^^y'~''^ , (14) 

'AxA 

assuming, as above, a polynomial bound in t for the commutators. 

In the following, the second and fourth momenta of du^Jk r ^^"^ 

their limits for R ^ oo will be calculated and discussed in terms 

of powers of fc, for A = p(/i), h{x) > 0. The one point result for 

the energy spectrum of plasma waves at zero momentum will use the 

following 

Lemma 2 If a sequence of positive even measures dvn{uj) satisfies 
Mm I dunuj"^^ = c^^ , A; = 0,1, 2 , c>0 , 



then the sequence converges in the sense of measures to 

l/2{5{uj-c) + 6{uj + c)) 
Iflim.nfdi'n = 1 and lim„ Jdf„u;^ = 0, then it converges to 6{uj). 

Proof: (tj^ — c^)'^ di'n{uj) are finite positive measures, their integrals 
converge to zero for n ^ oo and therefore, for any continuous bounded 
function /, J f{uj){uj'^ — c^)^ dun{uj) converge to zero. Any continuous, 
bounded, even function g can be approximated uniformly by 

<7(c) + /,(u;)(a;2-c2)2 
with /g continuous and bounded; this implies 

g{uj) dun -> g{c) 



and this is enough since the dun are even. The second statement 
follows similarly. D 
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3. Equations of motion and low k ex- 
pansions 

In the following, we analyze the energy spectrum associated to density 
waves in infinite ferniion systems. As a consequence of eq.Q, the 
calculation of the momenta of the corresponding measures dvpi}y\j{Lo) 
reduces to the application of the ACR relations to the time derivatives 
of p. By time translation invariance of $7, in order to obtain the 
momenta up to the fourth (the fifth for dp) , only three time derivatives 
must be computed. 

Assumptions A), B), C) will be needed only for the Wick poly- 
nomials of second degree, the only ones which appear in the time 
derivatives of p{x) = ifj* {x) iIj{x) (eqs.(14) and (19)-(24)). As usual, 
: M : will denote the Wick ordered polynomial M in the variables 
il:*{xi),ip{xi), i = 1...N \.e the polynomial with all the V'* on the 
left of all the ip and the sign of the corresponding permutation of the 
variables. For simplicity, we will omit the fermion field indexes. 

Considerable information about the low momentum behavour of 
charge density commutators can be obtained from general principles. 
The first observation is that each time derivative with respect to the 
free evolution of any second degree gauge invariant polynomial P{x) 
in the fermion fields explicitely introduces one more space derivative 
in the expression of P{x, t), and therefore one power of k in the Fourier 
transform of any correlation function of P. 

Lemma 3 For gauge invariant polynomials of second degree, P{x) = 
^On^kd^ '^*{x)d^'ip{x), d^""' denoting a product of n space deriva- 
tives, the commutator with the free Hamiltonian H^, in the limit 
R — > oo, is of the form ^^ diQi{x), with Qi in the same class. 

Proof: For the free evolution, from 

d/dt V'(x) = i lim [H%, '0(x)] = i/2m AV'(x) (15) 

R 

and the conjugate equation it follows immediately 

d/dtd^''^*{x)d^^^{x)= (16) 

= i/2m di{-did^''H*{x)d^^H{x) + d^''^*{x)did^^^{x)) D 

A second source of powers of k arises from symmetries which imply 
the vanishing of the integral of the commutators [A{x) , B] , for all local 
variables B, as operator valued distributions. In fact, by locality of 
the ACR relations such commutators have compact support, their 
Fourier transform is analytic in k and vanishes at A; = if the above 
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property holds. E.g., since the current ji generates space translations, 
the expectation value of current commutators on a space translation 
invariant state vanish at zero momentum: 



/ 



dxni[Ux),B])=0 (17) 



Prom parity invariance of 0, and definite parity of A, it also follows 
that 

{Aix)cr/deAiy)) 

is even and therefore its Fourier transform, if regular enough, is of 
order 2m + 2 in A; if it vanishes of order 2m. 

A basic result for the following is obtained from the commutator 
between the current 

ji{x) = i/2m {diip*{x) V'(x) - ^p*{x) dit{j{x)) 

and the interaction hamiltonian, which gives rise to the density of 
force 

Fi{x)=:p{x)jdydiVL{x-y)p{y): . (18) 

The commutator of Fi{x) with any local variable is integrable, also in 
the limit L ^ oo if diV{x) is integrable; it vanishes after integration 
in X since it is odd in the interchange of x and y (a consequence of the 
third law of Newton's). The same holds for the total force, for any 
number of fields. In the case of non integrable potentials, a subtraction 
is needed and a convenient relation is given by the following 

Lemma 4 For a space translation invariant state, if 

{[.{p{x)-{p))d,V{x-y)p{y):,B]) , (19) 

is integrable in the two variables, then its integral vanishes and 



dx lim / dy ([: (p(x) - (p)) diViix - y) p{y) : , B]) = . (20) 

In the absence of translation invariance, the same result holds with (p) 
replaced by {p{x)) and p{y) replaced by {p{y) — {p{y))) 

Proof: Using integrability, the l.h.s. of eq. (I20p gives the integral of 
()19p . by the Lebesgue dominated convergence theorem. On the other 
side, for all L, 

dy diViix - y) (p) = 



and therefore the l.h.s. of eq. (|2U|) is independent of the substitution 
of p{y) with p{y) — (p); in this form, the integral vanishes for all L by 
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antisimmetry in the exchange of x and y. The last statement follows 
from integrability and antisymmetry. D 

In the application below, the subtraction p{x) — (p) in Lemma 4 is a 
priori irrelevant for diV{x) integrable (and in fact also for didjV{x) 
integrable), but crucial for the Coulomb potential, where actually only 
the subtracted terms will survive in the analysis of the frequency spec- 
trum in the zero momentum limit. 

The equations of motion for the density are the following 

-pix,t) = -dmix) (21) 

^p{x,t) = l/m d, {dkSk,{x) + : p{x) di {Vl * p){x) :) (22) 

^P(.x,t) = +l/mdidk{Sl{x) +SUx)) + 

- l/m di (: dkjk{x) di {Vl * p){x) : + : p{x) di {Vl * dkjk{x) :) (23) 
with 

Sk^{x) = l/4m {dkiP* 5,V - dkd^r V + h.c.) , S^ = di Sik^ , (24) 

Siki being obtained as in Lemma 3 and 

SUx) = -{: jk{x) diVL * p{x) : +{i ^ k)) (25) 

In the presence of spin, exactly the same equations hold for p, jj 
and Sij replaced by the corresponding sum over the spin index. For N 
kinds of fermion fields, ipi, I = 1 . . . N with masses mi, the same equa- 
tions apply to the density operators p^'"', the corresponding currents 
j('' and stress tensors S^j^ , with mass m/ and the obvious substitution 

of diVL * p with Ylm ^i^L '™ * /O*-™"^- Since the potential is assumed to 
be regularized by the convolution with 7(x) G 5, and the correlation 
functions of p{x), ji{x), Sij{x) are assumed to be bounded after smear- 
ing, the space cutoff L can be removed in all correlation functions, for 
all potentials with integrable derivatives. 

Wick ordering in eqs. ([22]) . (p3]) . ([25]) can be omitted if the partial 
derivatives of the potential vanish at the origin up to the third order. 
Since this leads to simplifications in the calculations (only a few com- 
mutators are then needed, and one can forget about Fermi fields and 
their ordering) with no substantial consequences on the results [28], 
such a property will be assumed in the following for the regularization 
7 of the potentials (the derivatives of Vl vanishing at the origin up 
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to the third order if A7(0) = 0). Our calculations will only use the 
above time derivatives and the following basic commutators: 

[p{y),jj{z)] = {-i)/m d]{6{y - z)p{y)) (26) 

[p{y),S,u{z)] =-i{df{5{y-z)h{y)) + {i ^ k)) (27) 
[jj{y),Sik{z)] = -i/m{dl6{y - z)Sij(z) + {i ^ k)) + 

+ i/m d]{5{y - z) Sui{z)) - i/^m^ dfdld]{5{y - z) p{y)) (28) 



Eqs. (j26j) . (|27|) . (|28|) also hold for the sum over the spin index of the 
same operators. By eq.([9]), the required momenta are of the form 

= l/R''{-ir^\[p{h*f\n),d''+'/de+'p{h*fk,R,t)])\t=o . (29) 
For translation invariant states, 



/' 



^ a;" di^ito) = (27r)-'* / dq Tn{q) h{qY a^q " kY/R'' (30) 

with 

fn{k) ^ TTnix -y)^ i-ir+'T {[p{x), d-+^/de+^ p{y, t)] ) \t=o ■ 

(31) 
From the equations of motion and boundedness of the (equal time) 
correlation functions after smearing it follows immediately that Tn{k) 
is regular (C°° for potentials of fast decrease, continuous for potentials 
with diV{x) integrable), so that the limit i? — > cxo of eq. (|3U|) exists and 
gives Tn{k), apart from a constant times h{k)'^ , which is positive for 
small k. The corresponding normalized momenta of the frequency 
converge therefore to 

< a;"(A:) >= fn{k)/n{k) (32) 

For states invariant under lattice translations (omitting for sim- 
plicity the convolution with /i), 

l/R'' (-i)"+i / dx dy e'^'^y--^ aR{x) aR{y) T„(x, y) (33) 



with Tn{x, y) given by the r.h.s. of eq. (f3T|) . For potentials with diV{x) 
integrable, Tn{x,y) is integrable in y and its integral is periodic in x. 
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Therefore, the Umit i? ^ oo of eq. (|33p exists and gives, apart from an 
irrelevant constant, 

1/|A| j dx j dy e^'^^^"") T„(x, y) = fn{k, 0) , (34) 

with the integrand periodic in x, the notation referring to the Fourier 
series expansion in the variable x. Again, with the notation of eq.(|34 



<u;"(A;)>=T„(A:,0)/To(A:,0) (35) 

Eq.([26]) gives 

fo{k) = {p) k^/m (36) 



To(A:,0) = 1/|A| dx{p{x))k'/m={p)kym (37) 

Ja 

Eqs. (i2T]) . (f22l) immediately imply that T2{k) and T2{k,0) are in- 
finitesimal with respect to k'^, as a consequence of the presence of 
space derivatives (also following from Lemma 3) and of the arguments 
following eq. (|18|) . This also clearly applies for more than one kind of 
fermions, for the total density /9*°*(x) = "^Zi P i^)- Lemma 2 then 
gives a result similar to the Goldstone theorem for the spontaneous 
breaking of the Galilei group [27] , [9] : 



Proposition 1 For any number offermion fields and potentials V^^"^^ 
with integrable first derivatives, the energy spectral measures 
duptot(^ii^j^ ^, eqs. [T3\) ,p4^, normalized to total unit mass, converge as 
measures , for i? — > oo and A: — > 0, to S{u;). 

In the next Section the mean squared frequency is discussed for 
all potentials decaying at infinity faster than the Coulomb case; as 
a result. Proposition 1 will be extended to all such potentials, more 
precisely to all potentials with integrable second derivatives 

4. < uj'^{k) > and long range potentials 



The calculation of f2{k) and f2{k, 0) only requires eqs.jH]), ^, dM 
()28p . In fact, using time translation invariance of fi. 



T2{x,y) = -{-if{[d/dtp{x,t), d^/dt^p{y,t)])\t=o = 

-{-if/mdfd]{dl{[j,{x),SUy)]) + {[J,ix),p{y)dy{VL*p)iy)])) = 
= l/m2 d^Oyal {dl5{x - y) Sijiy) + . . . + 1/4 d^d^d]6ix - y) p{x)) + 

+l/rr? dfd] j dzdid,VL{y - z) {b{x - y) - 5{x - z)) {p{y)p{z)) + 
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+ l/m^ af A^(5(x - y){p{y) S^Vl * p {y))) (38) 



The dots refer to two terms obtained, as in eq. (|28p by permutations 
of indexes. In the presence of spin, each p and Sij operator must be 
summed over the spin index. The cutoff L can be removed if didjV{x) 
is integrable, as a consequence of boundedness (after smearing) of the 
two point correlation function of p. 

For translation invariant states, the last term in eg. (1380 vanishes by 
invariance under parity (or, see eq. (|22p . under space and time trans- 
lations), so that 

f2{k) = l/n? kikk (3 {Sik) k^-{Gik{k)-G,km) + l/Am^ {p) fc^ (39) 

with 

G,k{k)= T{dAV{x-y){p{x)p{y))) 

Eqs. (f30l) . (f36|) . ([39]) give, for the second momentum of the frequency 



< ^\k) > = ^^k,kk + A ^' - -^T^(G,k{k) - G,fc(0)) (40) 
In the rotation invariant case, 

{Sik{x)) = l/dm (d.rdjiP) = 2/d Ekra {p) 

Eq. ()40p has been derived in refs. [TT], [29] in the framework of a finite 
particle system in a large box (the use of the Coulomb potential in its 
derivation is not without problems, the plasma frequency being given 
by a singularity at the origin in Gij{k)). It holds, in our framework, 
for all potentials with didjV integrable. 

For states invariant under a lattice subgroup of translations, the 
last term in eq. ()38p vanishes by parity invariance after integration in 
X on a lattice cell; using periodicity in x and integration by parts, 

f2ik,0) = l/m^ kikki3{S,k)k^-iG',kik)-G'ikm) + l/^m' {p)k^ , 

(41) 
with p and Sik mean values over a lattice cell and 

Gl{k) ^ 1/lAI I dx j dye^'^^y-^^ dAV{x - y){p{x)p{y)) , 



so that, by eqs. (l35]) . (j37j) . the mean squared frequency is given again 
by eg. (Hoi) , with p and Sik substituted by their mean values and Gik 

by Gl. 

In all the cases, the low momentum behaviour of the mean square 
frequency depends on the decay property the potential at infinity. In 
fact, Gij{k) and G^j{k) are C" for potentials with \x\'^didjV integrable. 
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For n = 0, the (immediate) extension of eqs.(|40 p . (j4ip to A^ kinds of 
fermions proves Proposition 1 for potentials with integrable second 
derivatives. Actually, the result only depends on an additional power 
of the momentum appearing in all the commutators as a consequence 
of eq. (ll7p . For n > 4, in particular for sufficiently regular potentials 
decaying at least as \x\~'^~'^~'' , eq. (j4Up implies, using rotation invari- 
ance, 

< uj'^{k) >= {Q/dEkin + C) k^lm{p) + 0(k^) (42) 

The same holds for n >2 (in particular for regular potentials decaying 
at least as Ixl""^"^), with 0{k'^) replaced by o{k'^). In both cases C 
can be written, in dimension d = 3, 

C = 1/30 I d^x {p{<d)p{x)){x^/\ + 2xiXjdidj) V{x) 

For potentials with \x\^didjV{x) integrable, < /3 < 2, the last term 
in eq. ()40p can dominate at low momentum; it can be estimated as: 

\G^J{k)-Gi,{0)\ <2\k\f' fdx \x\f^\Gi,{x)\ , 

so that 

< LO^{k) > < const \kf . (43) 

If the potential term dominates, the low momentum behaviour of 
< uP'ik) > is clearly independent from the addition of any potential W 
with x^ didjW integrable; for potentials V{x) ~ |a:|~°, d — 2<a<d, 
one obtains 

<UJ^{k) >~ \kf 

with /? = a — d + 2. The same results apply to states invariant under 
lattice translations. 

In the following Sections the fourth momentum of the frequency 
is discussed, in the zero momentum limit for the Coulomb interaction 
and to the order k'^ for short range potentials. 

5. The plasmon spectrum in Jellium 

We will derive in this Section the one point result, eq.([T]), for the en- 
ergy spectrum of density waves in the limit A; ^ for the jellium 
model, in three space dimensions; the result will follow from the ap- 
plication of Lemma 4 to the analysis of the fourth momentum of the 
frequency and of Lemma 2 to the momenta up to the fourth. 

As discussed above, in the case of a Coulomb potential an infrared 
cutoff is necessary in the equations of motion. Its removal in the first 
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four momenta of the frequency only requires that the two and three 
point correlation functions of p{x) — {p{x)) decay (respectively) as 
\xi — X2\~^~'^ and |xi — X2\~^~'^ \x2 — xsj^-^^^. 

Eqs. (l40p . (l4ip for the second momentum of the frequency holds for 
potentials Vl with ultraviolet and infrared regularization, eqs.([6]),([7]) 
The first two terms in eq. (|4U|) . are of order k"^ and /c^; for the last 
term, from integrability in y of eq. 

djV{x-y){{p{x)p{y))-{p{x)){piy))) 

it follows that, after such a subtraction, the limit for L ^ oo vanishes 
of order k. We are left therefore with 

TididjVLix-y){pix)){p{y))) 

and therefore, in the space translation invariant case, 

< uj'^(k) >= lim A;2 VL{k){p)/m + 0{k) -^ e^{p)/m (44) 

for /c — > 0, using 7(0) = 1. 0{k) becomes 0{k'^) if the truncated 
correlation function of p decays as Ixl"^"*". 

For lattice translation invariant states, expanding the periodic 
function {p{x)) in Fourier series, the corresponding term reads 

kk- v^ 

— hrn? /. I iPn)\^ {{ki + ni){kj + Uj) Viik + n) - UiUj Viin)) , 
m (p) fe^ ^-^ 

n ranging over the reciprocal lattice of the translation lattice. For 
L — > cxD, all the terms with n ^ are of order k and the same holds 
for their sum (as a consequence of the regularization of the potential) . 
Therefore < co'^^k) > converges, for A; ^ 0, to e'^{p)/m (in fact p = po), 
i.e. to the square of the plasma frequency associated to the mean 
density. 

The absence of two powers of k in eq. ()44p with respect to eq. 
()42p depends on the failure, for the Coulomb potential, of both the 
mechanisms mentioned in Sect. 3, i.e., in the subtraction needed for 
the validity of Newton's third law and in the non-integrability in x of 
the commutator 

-i / dy {[ji{x),dkV{y)p{y)]) = -1/m dkdiV{x) {p{x)) , 

which in fact gives, when summed over equal indexes, e^(p) (Gauss' 
law). The result for the mean squared frequency also follows from an 
analysis of the removal of the infrared cutoff in the above commutator: 



-\/m I dx didi Vl{x) p{x) = e /m j dx {S{x) — aL{x))p{x) 
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with aL{x) of support near \x\ = L and J uiix) dx = 1, so that, in all 
correlation functions, for L -^ co, J dx (Jl{x)p{x) -^ (p) 

The fourth momentum of the frequency is obtained, eqs. (j3ip . (j32p . 
from the commutator between the r.h.s. of eqs. (l22|) . ()23p . The com- 
plete commutator has been calculated in ref. [28j. In order to discuss 
the zero momentum limit, by eqs. ([S^ . ([5B]) . (f57|) . we need only terms 
up to the second order in k, i.e. the commutators between the last 
terms in equations (j22p . (j23p . We need therefore to discuss the limit 
L ^ oo and then A; — > of 

( [p{x) diVL * p{x) , dkjk{y) dj Vl * p{y) + p{y) dj Vl * dkjkiv)] ) (45) 

By applying eq. (|26p . four terms appear, each involving the three point 
function of the density. In the translation invariant case, y can be fixed 
and (p) subtracted in the convolutions. If we also subtract to the l.h.s. 
of the commutator a term (p) diVL*p, such terms are integrals, in two 
variables, of functions bounded by integrable functions uniformly in L, 
as a consequence of the decay assumptions on the correlation functions 
of the density. Their limit L — ;• oo vanishes by Lemma 4 and therefore 
we can substitute {p) to p{x) in the l.h.s. of the commutator. The 
same argument then applies, by antisymmetry of diV{z) to the r.h.s., 
for the subtraction p{y) — > p(y) — (p) We are therefore left with 

df d]{ [ {p) diVL * p{x) , {p) d, Vl * duUv) ] ) 
and eq. (p6|) gives 

(p)2 f dw dz AVLix-w) AVL{y-z){-i)/mdld'^{5{w-z)p{w)) = 

= -i/m {pf JdzdkA Vl{x -z)dkA V^y - z) , 

so that, by eq.([HT]). 

J'Tiik) ~ l/rrV^ {pfik^f Vl{k) ^ l/rr? {pf e^ k^^k) 

for L ^ oo. This immediately implies 

< uj\k) > = fi{k)/n{k) -^ e^ {pf /m^ (46) 

for A; — > 0, independently of the regularization 7, with 7(0) = 1. 
By Lemma 2, eqs. (jM]), (fl6l) imply the one point result, eq.(l), for 
the energy spectrum of density waves in the zero momentum limit. 
The result is independent from the ultraviolet regularization of the 
potential and clearly holds for all potentials with k'^V{k) -^ e^ for 
k^O. 
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For states invariant under lattice translations, eqs. (j34p .(|35 p apply 
and we must consider the mean in x and integral in y of the commu- 
tator (j45p . Using the cluster properties and applying Lemma 4, the 
commutator vanishes, for L — > oo, at zero momenta, after the sub- 
tractions p ^ p — {p{x)) in the l.h.s; the remaining terms vanish, as 
above, after the same subtraction in the r.h.s., so that we end with 
the expectation value of the commutator between 

{p{x)) {d^VL * p){x) + p{x) {{d^VL * p){x)) 

and 

dkjk{y) {{dj Vl * p){y)) + {p{y)) dj Vl * dkjkiv) 

Using integration by parts in eq. ()34p and discarding terms of order 
> k'^ (uniformly in L), eq. (|26p gives 



T4{k,0) ~liml/m |A| dx dydw dzdidkVL{x — w)djdkVL{y — z) 

kikj e^^(s^-^') p^p^ {py {S{w-z)-5{x-z))+p^ {5{x-y)-6{y-w))) (47) 

with Pa = (p(o)) In the translation invariant case, only the first term 
appeared in the r.h.s., the others corresponding to an irrelevant sub- 
traction of a constant in the convolutions. By the previous result on 
the second momentum of the frequency, the one point result for the 
plasma spectrum is equivalent to the cancellation of all the terms dif- 
ferent from the mean in the Fourier expansion of the expectation value 
of the density. The expansion gives, with pn the Fourier coefficients 
oi Px, 

hm ^ pni pn2 Pn-i {didkVL{-k - ni) - didkVL{n2)) 

ni+n2+n3=0 

{djdkVL{k - ns) - djdkVLins)) ; 

by the regularity of V{k) for A; 7^ and the vanishing of djdkVL{0)) , 
all the terms are of order k uniformly in L, except those with ns = 0; 
in this case, since ni = —n2, the limit is non vanishing only for ni = 
^2 = ^3 = 0. By the fast decrease of V{k) the same holds for the sum. 
The result is therefore, at the order k'^ 

nik, 0) ~ {pf/m^ {k^f V{-k) V{k) ~ e^ {pf k^ jm^ ; (48) 

as before, eq. (l48p holds independently of the ultraviolet regularization 
7, only requiring k'^V{k) — > e^ for k ^ 0, and implies eq. (j46p . with 
(p) replaced by the mean density {p). In the presence of spin, the 
same derivation and results apply to the charge and current operators 
summed over the spin index. 
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6. < LJ^{k) > at the order /c^ 

We calculate the fourth momentum of the frequency at the order k'^ , 
for translation and rotation invariant states, from eqs. (|31|) . (j32p . From 
Sect. 3 it follows that, for potentials of fast decrease, Ti{k) is regular 
(C°°) and at least of order k'^. In ref.|2H] Ti{k) has been obtained, 
at the order A;^, in terms of two point functions {pSik), {jijk) and 
of the three point function of p. Using identities which follow from 
the equation of motion, we shall express T^ik) at the order /c^, and 
therefore < Lu'^{k) > at the order k'^, purely in terms of the three point 
function of the density, for one component models with short range 
potentials. By eq. ([3T]) . 



T^ix -y) = i-i) {[d^/dt^pix, t), d'/dt^ piy, t)] ) (49) 

at t = 0. In the commutator of the r.h.s. of eqs. ()22p . (l23p we drop 
Sfj^{x), which produces terms of order k^ (also as a consequence of the 
cancellations discussed in Sect. 3), use the identity (for A and B of the 
same definite parity) 

{[A{x),B{y)]) = {[B{y),A{x)]) = {[B{x),A{ym 

and eq. ()25p . We obtain, using translation invariance, 

ini^ Ti{x -y) = - 2didk dj dm{[Skiix) , jj{y) dmV * p{y)]) 

- 4 di dj d„,{[p{x) d,{V * p){x) , j,{y) draV * p{y)]) (50) 

+ did,{[p{x)di{v*p){x), {dkjk{y)d,{v*p){y)+p{x)dj{v*dkjk{y))) 

where d denotes the derivation with respect to x. The second and 
third commutators only require eq. ()26p . so that the result involves 
purely the three point function of p. The first commutator requires 
eqs. (f28]l and (f271) : dropping from eq. (f28]) the last two terms, of higher 
order in k, omitting the four derivatives and integrating in y we obtain, 
for the coefficient of kikjk^km, 



4i / dz dkdmV{x - z) 1/m {Sij{x)p{z) - ji{x) jj{z)) 

Rotation invariance, symmetry under the permutations i *^ j and 
i '^ k and summation over equal indices give, for the coefficient of A:^, 



4V5 / dz V{x - z){l/m {didjSijix) p{z)) + {diji{x) djj,{z))) 

The identity 

= d/dt{p{x,t)p{z,t))\t=o = {diji{x)djjj{z)) + 
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+ 1/m {di djSij{x) + di{p{x) diV * p{x)) p{z)) 
then gives, for the first term in T4, 



Tl'\k) = -4/5 kym' / dydzd^V{y)d^V{z){p{Q)p{y)p{z)) (51) 



The second and third commutators in eq. (j50p are immediately cal- 
culated from eq. ()26p : in both commutators, a term with four space 
derivatives appears, of the same form as in eq. (|5ip . with coefficients, 
respectively, 4/3 and 1/3. The remaing terms only contain three or 
two space derivatives, and a Taylor expansion in k of the correspond- 
ing expressions is necessary; the resulting contributions to T4^{z) are 
therefore integrals of the three point function of the density with first 
and second order polynomials in z. Summing all the terms, we obtain, 
to the fourth order in k. 



T4k) = (-4/5 + 5/3) kym-' I dy dz diV{y) d,V{z) (p(0) p{y) p{z)) + 
+ kikjkikm/m^ dy dz Zijimiy,z) {p{0) p{y) p{z)) , (52) 



with 

Ziji.m{y,z) = 6diV(y)zmdjdiV{z)+ 

+1/2 didkViy) zi zm djdk{2V{z) - V{z - y)) . 
By rotation invariance, the result can also be written as 

n{k) = k^rr? I dy dz W{y, z) (/^(O) p{y) p{z)) (53) 

with 

W{y, z) = 13/15 diV{y) diV{z) + 6/15 d,V{y){ziA + 2z-d d,)V{z)+ 

+ 1/15 didkV{y) djdkV{z) {6ij y ■ z + yiZj + yj Zi)) (54) 

and eq.(l2|) follows. 
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